Non-Koopmans Corrections in Density-functional Theory: Self-interaction Revisited 
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In effective single-electron theories, self-interaction manifests itself through the unphysical depen- 
dence of the energy of an electronic state as a function of its occupation, which results in important 
deviations from the ideal Koopmans trend and strongly affects the accuracy of electronic-structure 
predictions. Here, we study the non-Koopmans behavior of local and semilocal density-functional 
theory (DFT) total energy methods as a means to quantify and to correct self-interaction errors. We 
introduce a non-Koopmans self-interaction correction that generalizes the Perdew-Zunger scheme, 
and demonstrate its considerably improved performance in correcting the deficiencies of DFT ap- 
proximations for self-interaction problems of fundamental and practical relevance. 
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Because density-functional theory (DFT) accounts for 
correlated electron interactions via an explicit functional 
E xc of the total density p, it provides a computational 
scheme that is far more predictive and practical than 
the Hartree-Fock (HF) method for determining the Born- 
Oppenheimer energy of quantum systems [lj . In spite of 
the proven accuracy of local and semilocal DFT total en- 
ergy calculations, it has long been recognized that the 
eigenvalues e, a of the Kohn-Sham Hamiltonian, which 
are defined in the strictest sense as the Lagrange multi- 
pliers associated to the orthonormality constraint in the 
self-consistent energy minimization, have limited phys- 
ical relevance 2]. Indeed, due to the fact that Koop- 
mans' theorem does not hold for approximate DFT func- 
tionals, the Kohn-Sham single-electron energies are not 
related to any physical ionization process, at variance 
with their HF counterparts that can be directly identified 
as the opposite electron removal energies -AE^J of the 
unrelaxed electronic system (i.e., without self-consistent 
modification of the electronic wavefunctions) . Following 
the seminal work of Perdew and Zunger Q , one can ac- 
curately quantify Kohn-Sham deviations from the ideal 
Koopmans behavior by defining non-Koopmans energy 
contributions n iQ = a a + A£j|°' (this heuristic definition 
will be made explicit later in the course of our discussion) . 

As pointed out by Perdew and Zunger, for any self- 
interaction-free single-electron theory, the Koopmans 
condition 



U ia = 



(1) 



is rigorously satisfied (in particular, the HF functional 
verifies 11™ =0 precisely because of the cancelation of 
the Hartree and exchange self-energy terms) Q. This 
non-self-interaction criterion establishes a central and ex- 
act correspondence between the non-Koopmans behavior 
of a given energy functional and self-interaction errors, 



which favor electronic delocalization and considerably af- 
fect the accuracy of local and semilocal Kohn-Sham cal- 
culations. Indeed, self-interaction is responsible for well- 
known quantitative and qualitative failures of conven- 
tional DFT functionals in describing important phenom- 
ena, such as electron transfer [3] , electronic transport [|| , 
and electrical polarization in extended systems [6(. 
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FIG. 1: HF and LSD Born-Oppenheimer hydrogen total en- 
ergies E and single-electron energies e with and without elec- 
tronic relaxation as the function of the fractional occupation 
of the Is orbital. 



In this work, we examine the cancelation of the non- 
Koopmans terms Hi a that results from the Perdew- 
Zunger self-interaction correction (PZ-SIC) and show 
that the non-self-interaction condition (Eq. [I} is satis- 
fied to second order in the single-electron densities. In 
order to achieve a higher level of accuracy, we introduce 
an alternative non-Koopmans self-interaction correction 
(NK-SIC). We validate our method by studying the dis- 
sociation of H^ and H2 molecules, and compare its per- 
formance to that of the original PZ-SIC approach in de- 
scribing the longitudinal electrical response of dimerized 
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hydrogen chains, a complex self-interaction problem that 
has received much attention for its fundamental relevance 
in molecular optoelectronics [f| 0, S, 0, E3| ■ 

As a prelude to studying self-interaction in many- 
electron systems, we analyze the non-Koopmans behavior 
of the local spin density (LSD) functional in the simplest 
case of the H — ► H + ionization (Fig. [I]). The dependen- 
cies of the relaxed energy £; LSD and unrelaxed energy 
_gLSD,(o) as a f unc tion of the fractional occupation / of 
the LSD Is state are depicted and compared to HF in 
Fig. [U (in these calculations, we use the standard gen- 
eralization of the HF theory to fractional occupations). 
The derivatives of the LSD energies, which correspond to 
the relaxed and unrelaxed single-electron energies e LSD 
and e LSD '(°) (Janak's theorem), are also shown. Since, 
for this one-electron system, the relaxation contribution 
e lsd = ^lsd _ £Lsd,(o) vanis h es at / = and / = 1, 
the ionization energy AE LSD = E LSD (0) - £ LSD (1) is 
exactly equal to its unrelaxed counterpart Ai? LSD '(°). 
Note however that the small relaxation energy vjLSD re _ 
mains strictly negative in the transition state region at 
variance with HF, for which no relaxation occurs. Al- 
though LSD predicts the energy of the hydrogen atom to 
be —0.94 Ry in agreement with the exact HF result of —1 
Ry, the LSD energy departs significantly from the linear 
Koopmans behavior due to unphysical self-interaction: 
the slopes e LSD and e LSD <(°) are largely underestimated 
close to / = 0, reflecting the predominance of the neg- 
ative 0(/ 4 ^ 3 ) exchange-correlation contribution, while 
they are overestimated around / = 1 due to the pos- 
itive 0(/ 2 ) Hartree term. Despite these notable non- 
Koopmans deviations, we observe that the LSD energies 
approach the expected Koopmans behavior at / = 



tional and orbital-dependent Hamiltonian are defined as 



(i.e., e 



LSD,(0) 



,LSD 



-1 Ry). As an im- 



portant consequence, the non-Koopmans contribution to 
the unrelaxed ionization energy must be defined by ref- 
erence to the single- electron energy of the empty state: 
n = e(°)(0) + AE(°\ Extending these observations to 
many-electron systems with fractional occupations and 
making use of Slater's theorem [3|, the non-Koopmans 
self-interaction energies can be rewritten as 



Hi. 



rfta 

dX 



AE, 



(0) 



ii 



(42(0) 



(0) 



(A) 



(2) 



where e\ a ' (A) denotes the effective energy of the orbital 
ipi a with occupation A in the unrelaxed transition state. 
Note that in the absence of self-interaction, the unre- 
laxed single-electron energies satisfy de^ a '(X)/dX = (by 
definition) and Eq. [I] is verified. 

We are now in a position to evaluate the non- 
Koopmans energy cancelation resulting from PZ-SIC, 
which consists of subtracting single-electron Hartree and 
exchange-correlation contributions to the Kohn-Sham 
functional E KS . Explicitly, the Perdew-Zunger func- 



E 



PZ 



E KS -J2 ia E H [p ia ]+E xc {p ia ] 



h ia = h - VHlp ia ] - v KCt0t [p ia ], 



(3) 
(4) 

where pi a — fi a \ipi a \ 2 denotes single-electron densities, 
vn = 5 En/ 5 p is the electrostatic Hartree potential, and 
w XC) a = SE xc /5p a stands for the exchange-correction po- 
tential. Rewriting the non-Koopmans energies (Eq. ^ 
in terms of the corrected Hamiltonian (Eq. [4| and ex- 
panding each energy contribution in terms of the vari- 
ables pi a at the self-interaction-free empty-state density 
Tha = (Pa — Pia, Pa) (with a opposite to a), we obtain: 



n 



PZ 



10 

£ 



dX{^ ia \(h^(0)-h^(X)) I Via) 

J f^a aa [p~^}PjaPi a Pi a dridr 2 dr 3 , (5) 



where /i™^ [p] denotes the nth-order functional derivative 
of E xc [p\. As a result, PZ-SIC cancels self-interaction 
contributions up to the second order in the single-electron 
densities, which considerably ameliorates the precision 
of calculated ionization energies and related properties 
0, [ll[ but tends to sensibly overestimate self-interaction 
errors for many-electron systems 12|, |l3|, UA | • 

A more accurate self-interaction cancelation can be 
obtained by modifying the expression of the corrective 
terms in Eq. [3J In fact, due to the direct correspon- 
dence between non-Koopmans contributions and self- 
interaction errors, it is legitimate to identify the or- 
bital self-interaction correction as the single-electron non- 
Koopmans energy n^ s . The resulting NK-SIC energy 
reads E NK = E KS + J2 ia ^fa an d variation with re- 
spect to the single-electron density pi a yields the orbital- 
dependent Hamiltonian 



hl K = h 



KS 



[p ia\ T" ^xc.a 
V -f- 1?H [Pia] ~t~ ^xc,a[pm] ~r~ ^xc, 



(6) 



where v is the external potential and tv^cia 

is defined as 
w xc , ia (r) = Ej^fa^xc.aC^sr) -u xc , a (M;r) 

+ / dr>f£i p (lpJe}-ry)p jP (T>). (7) 

Substituting Eqs. [6] and [7] in the expression of the non- 
Koopmans energies (Eq. [3]), we obtain the following es- 
timate for the self-interaction contributions: 



HL K x ^ 



^_ , /ic,L/3/3([Pi/3]; r l J r 2 I I'3,r4) 

' Pjf3 {^i)pjj3 (f3)Pia fa )pia (r 1 )dY 1 dY 2 dr 3 dr 4 



(8) 



(note that only the term w xc .i a contributes to self- 
interaction). Consequently, our modification of the PZ- 
SIC functional brings about an additional factor of preci- 
sion in the cancelation of spurious electronic interactions. 
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FIG. 2: Exact, NK-SIC(LSD), and LSD potential energy 
curves of as a function of the H-H bond length. 



The improved accuracy of the NK-SIC scheme in com- 
parison to PZ-SIC will be shown below. 

At this point, one may observe that an exact cancela- 
tion of the non-Koopmans self-interaction contributions 
is obtained by defining the orbital-dependent Hamil- 
tonian as h ia = -5V 2 + v + v- H [p~] + w xc ,q[p^] = 
— 7;V 2 + Vi a [pi^). Nevertheless, such a scheme does not 
derive from the minimization of any physical energy, 
as this would lead to the following variational inconsis- 
tency: » £g = 5^ ± gsia = d'E The NK _SIC 
Hamiltonian (Eq. [5]) , which includes the additional term 
0(J2jf3^ia \\pjp\\) for variationality, provides a 
close physical approximation to this exact correction. 

The implementation of the NK-SIC method requires 

computing the kernel of the local or scmilocal cxchangc- 
(2) 

correlation energy /xc ■ This calculation is performed 
using the scheme proposed by Dal Corso and de Giron- 
coli in the context of phonon-dispersion computations 
[lit . Moreover, the determination of the NK-SIC terms 
necessitates calculating the Hartree self-interaction en- 
ergy En[pi a ] and potential i>H[pia], which entails treat- 
ing a system with a net electrical charge. To elimi- 
nate periodic-image errors in the plane-wave evaluation 
of En[pi a ] and we employ countercharge correc- 

tion methods [161 ]. Since the determination of /x C and 
the periodic-image corrections are relatively inexpensive, 
the computational burden of NK-SIC computations is 
comparable to that of conventional PZ-SIC calculations. 

We first study the dissociation of the Hj" molecular 
ion. The LSD and NK-SIC (LSD) energies are compared 
to the exact Born-Oppenhcimcr potential energy curve 
in Fig. [2j We observe that the LSD energies of this one- 
electron system are largely overestimated by as much as 
0.1 Ry in the vicinity of the equilibrium bond length due 
to self-interaction errors, at variance with the NK-SIC 
energies, which are in close concordance with the exact 



FIG. 3: UHF, NK-SIC(LSD), and LSD potential energy 
curves of H2 in the absence of electronic correlation as a func- 
tion of the H-H bond length. 

results. Moreover, in the infinite separation limit, the 
NK-SIC electronic ground state is correctly predicted to 
be degenerate (with an energy of —1 Ry) with respect 
to the effective occupations n\ and ria = 1 — ni of the 
two hydrogen sites, at variance with LSD, which yields a 
nondegenerate and delocalized split-electron state {n\ = 
712 = 1/2) of lower energy. 

As a second self-interaction paradigm, we examine the 
dissociation of H2 in the uncorrelated limit. The de- 
pendencies of the unrestricted Hartree- Fock (UHF), un- 
correlated NK-SIC (LSD), and uncorrelated LSD Born- 
Oppenheimer energies as a function of the intramolecu- 
lar distance are depicted in Fig. [3J Similarly to the H^ 
case, we observe that LSD overestimates the total en- 
ergy of the system, in contrast to NK-SIC, which brings 
the potential energy curve in close accordance with the 
sclf-interaction-free result. Furthermore, the analysis of 
the spin occupations n; Q of the hydrogen sites at large 
atomic separation reveals that both the NK-SIC and LSD 
methods predict a two-fold degenerate ground state, in 
which each electron is localized at one site (n^ = n2 [ = 1 
or nil = n 2| = !)■ However, this transition occurs for 
distinct reasons in the two cases: while electronic local- 
ization is driven by the reduction of the interelectronic 
Hartree energy for NK-SIC, it is the enhancement of un- 
physical self-exchange that leads to localization for LSD. 
This fact explains salient differences in the asymptotic 
behaviors of the LSD and NK-SIC energies. 

Having validated the NK-SIC method for archetypical 
one- and two-electron problems (for which PZ-SIC would 
also be exact), we turn to the difficult case of the elec- 
trical response of extended hydrogen chains, which are 
frequently employed as model systems in molecular op- 
toelectronics and represent a critical test in assessing the 
predictive performance of electronic-structure methods 
[17| . We consider a linear chain geometry characterized 
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TABLE I: Longitudinal polarizability ajv (in a.u.) of the 
hydrogen chain as a function of the number of dimers N. 
The mean absolute deviations Aa from CCSD(T) results are 
also reported. 



method 


OL2 


Q3 


Q'4 


Q6 


Aa 


CCSD(T) a 


29.0 


50.9 


74.4 


124.0 




MP4 6 


29.5 


51.6 


75.9 


126.9 


1.4 


NK-SIC(LSD) 


30.9 


49.0 


71.3 


122.3 


2.1 


HF" 


32.2 


56.6 


83.0 


138.6 


8.0 


PZ-SIC(LSD) a 


33.0 


59.7 


89.1 


152.0 


13.9 


LSD C 


37.6 


72.9 


115.3 


212.1 


39.9 



Ref. [13, 6 Ref. H, c Ref. [|. 



tions in comparison to PZ-SIC. In order to validate the 
NK-SIC method, we have studied the dissociation of H^ 
and H2, finding excellent agreement with exact and self- 
interaction-free calculations. Finally, we have demon- 
strated the much improved performance of the NK-SIC 
scheme in describing the electrical response of hydrogen 
chains, which represents a stringent and technologically 
relevant benchmark for electronic-structure methods. 

The computations in this work have been performed 
using the Quantum-Espresso package The authors 
acknowledge support from the Grant in Aid for Research, 
Artistry and Scholarship of the Graduate School of the 
University of Minnesota. Helpful suggestions and com- 
ments from Stefano de Gironcoli, Eric Cances, Nicolas 
Poilvert, and Andrea Floris are gratefully acknowledged. 



by alternating H-H distances of 2 and 3 bohr, as reported 
in the literature [2(j. The LSD and NK-SIC longitudinal 
polarizabilities a at of the hydrogen chain as a function of 
the number of H2 units are compared with HF, fourth- 
order M0ller-Plesset perturbation theory (MP4), coupled 
cluster [CCSD(T)], and PZ-SIC predictions in Table H 
We observe that the computed polarizabilities follow an 
increasing and progressively linear trend as a function of 
the number of hydrogen dimers. In spite of these sim- 
ilarities, we note that LSD dramatically overestimates 
the electrical response of the system. The Perdew-Zunger 
correction sensibly improves the accuracy of LSD results, 
reducing the mean absolute error from Acv LSD = 39.9 to 
Aa pz = 13.9 a.u. relative to CCSD(T) polarizabilities 
(which we choose as our standard of accuracy). Despite 
this notable amelioration, the performance of the PZ-SIC 
scheme remains inferior to that of self- interaction-free HF 
(Aa HF = 8.0 a.u.). This should be contrasted to NK-SIC 
predictions, which are found to be much more accurate 
than their HF counterparts with a precision comparable 
to the mean absolute difference between post-Hartree- 
Fock CCSD(T) and MP4 calculations (Aa NK = 2.1 a.u.). 
These findings confirm the self-interaction origin of the 
overestimated polarizabilities within LSD and provide a 
clear illustration of the superior predictive ability of the 
NK-SIC method for complex self-interaction problems. 

As a final note, it is important to mention that several 
successful refinements of the PZ-SIC method have been 
proposed in the literature [l2|, [HI- One limitation of 
such approaches is that they require to scale down the 
Hartree self-energy although this corrective term is in 
principle exact (as a result of the linearity of «h [/>])• The 
NK-SIC method represents a parameter-free alternative 
to semi-empirical downscaling procedures. 

In conclusion, we have revisited self-interaction in light 
of the Koopmans condition (Eq. [1]) and introduced a 
non-Koopmans correction (NK-SIC), which provides a 
more accurate cancelation of self-interaction contribu- 
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